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We propose a new approach to study quantum phase transitions in low-dimensional fermionic or 
spin models that go from uniform to spatially inhomogeneous phases such as dimerized, trimerized, 
or incommensurate phases. It is based on studying the length dependence of the von Neumann 
entropy and its corresponding Fourier spectrum for finite segments in the ground state of finite 
chains. Peaks at a nonzero wave vector are indicators of oscillatory behavior in decaying correlation 
functions. 
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Recently, it has been shown []J that, quite generally, 
discontinuities in some measure of entanglement between 
different parts of a ground-state quantum system can be 
used to study quantum phase transitions (QPTs). This 
measure can be the concurrence l2j , which can be used to 
study spin- 1/2 models 0, HH0,H 0,0, the entropy of a 
single site [llillllll2|> the entropy of a pair of neighboring 
sites |l3 T o r, more generally, the entropy of a block of 
length I [HEi. 

In a parallel development, it has been pointed out that 
critical and noncritical systems behave differently 0,0] 
when the length dependence of the entropy of a finite seg- 
ment of a quantum system is studied. The von Neumann 
entropy of a subsystem of length I (measured here and 
subsequently in units of the lattice constant a) saturates 
at a finite value when the system is noncritical, i.e., when 
the spectrum is gapped, while it increases logarithmically 
for critical, gapless systems. An analytic expression has 
been derived for models that map to a conformal field 
theory 0] , and this form has been shown to be satisfied 
by critical spin models. The entropy for a subsystem of 
len gth I in a finite system of length N has been shown to 
be 13 
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where c is the central charge and g is a constant shift due 
to the open boundary which depends on the ground-state 
degeneracy (lij . 

The aim of this paper is to show that the length de- 
pendence of the von Neumann entropy of a subsystem 
can, in fact, display a much richer structure than dis- 
cussed until now, and that its analysis allows a better 
characterization of the QPT, thereby providing a new 
diagnostic tool to study QPTs. Moreover, this method 
seems to be appropriate to study cases when no true 
phase transition takes place, i.e., when only the char- 
acter of the decaying correlation function changes. Our 



method is especially convenient when the density-matrix 
renormalization-group (DMRG) algorithm |2fjj is used, 
in which the density matrix of blocks of different lengths 
are generated in the course of the procedure so that the 
von Neumann entropy can be easily calculated. 

The first, simplest case which we will consider is the 
Majumdar-Ghosh (MG) model [2l|, which corresponds 
to the frustrated spin- 1/2 Heisenberg chain, described 
by the Hamiltonian 



H = J2 [ J (Si • S i+1 ) + J' (St ■ S l+2 )] 



(2) 



at the particular parameter value J' = 0.5 J. Since the 
ground state (of a finite chain) consists of independent 
nearest-neighbor singlets between odd and even sites, the 
entropy of a segment of length I oscillates between In 2 
when the block has an unpaired spin at its end and 
when all spins in the block are paired into singlets [2l| . 
We will also consider the exactly solvable Takhtajan- 
Babujian (TB) 22] and Uimin-Lai-Sutherland (ULS) 
points of the spin-one bilinear-biquadratic model 24], 



H = Y J [vos6{Si ■ S i+1 ) + sin6(Si ■ S l+1 ) 2 



(3) 



which correspond to the parameter values 9 = — 7r/4 and 
7r/4, respectively. As shown in Fig. ^ a periodic oscil- 
lation is superimposed onto a curve that is described by 
the analytic form given by Eq. in both cases. At the 
TB point, the period of oscillation is two lattice sites, 
while at the ULS point it is three lattice sites. When 
the length I is taken to be a multiple of two for the TB 
point or a multiple of three for the ULS, the entropy s(l) 
can be well-fitted using Eq. (JIJ with c approaching the 
known values, c = 3/2 J2f| and c = 2 |2g, respectively, 
in the limit of large N. 

In order to analyze the oscillatory nature of the finite 
subsystem entropy s(Z), it is useful to consider its Fourier 
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FIG. 1: Length dependence of the von Neumann entropy of 
segments of length / of a finite chain with N = 60 sites for (a) 
the Takhtajan-Babujian and (b) the Uimin-Lai-Sutherland 
models. The solid lines are our fit using Eq. Q taking every 
second and third data points. 
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with s(0) = s(N) — where q — 2-KnjN and n — 
—N/2,...,N/2. Since s(l) satisfies the relation s(l) — 
s(N—l), its Fourier components are all real and symmet- 
ric, s(q) = s(— q)] therefore, only the < q < tt region 
will be shown. Except for the large positive s(q = 0) com- 
ponent that grows with increasing chain length, the other 
components are all negative. They are shown for the two 
cases discussed above in Fig. As can be seen, apart 
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FIG. 2: Fourier spectrum s(q) (scaled by the system size N) 
of the length-dependent von Neumann entropy of finite chains 
of length N — 30, 60 and 90 for (a) the Takhtajan-Babujian 
and (b) the Uimin-Lai-Sutherland models. 

from the q — point, the Fourier spectrum exhibits (neg- 
ative) peaks at q = 7r and q = 2n/3, respectively. This 



is related to the fact that the TB model has two soft 
modes, at q = and 7r, while the ULS model has three, 
at q = and ±27r/3. Although finite-size extrapolation 
shows that these components vanish in the N — * oo limit, 
these peaks in the Fourier spectrum are nevertheless in- 
dications that the decay of correlation functions is not 
simply algebraic in these critical models, but that the 
decaying function is multiplied by an oscillatory factor. 
When the same calculation is performed for 8 in the range 
— 3tt/4 < 8 < — 7r/4, where the system is gapped and 
dimerized, the peak at q = 7r remains finite as N — > oo. 
On the other hand, in the whole interval ir/4 < 8 < ir/2, 
where the system is gapless and the excitation spectrum 
is similar to that at the ULS point, the entropies for block 
sizes that are multiples of three can be well-fitted with 
the form given in Eq. Q with c = 2, and for finite chains, 
a peak appears in s(q) at q = 2tt/3, in agreement with 
Refs. |3 an d Hil Thus, peaks in the Fourier spectrum 
of the length-dependent block entropy can provide useful 
information about QPTs, even when they scale to zero 
in the thermodynamic limit. 

We demonstrate this procedure on the example of the 
spin-one bilinear-biquadratic model near the VBS point 
|27j |. corresponding to #vbs = arctanl/3 ~ 0.10247T. It 
is known |28l | that this point is a disorder point, where 
incommensurate oscillations appear in the decaying cor- 
relation function; however, the shift of the minimum 
of the static structure factor appears only at a larger 
9j, = 0.1387T, the so-called Lifshitz point. In earlier work 
|l3j, some of us showed that s(N/2) has an extremum 
as a function of 8 at 6*vbs- Here we show that this ex- 
tremum is the indication that, in fact, #vbs is a dividing 
point which separates regions with a different behavior 
of s(l) and s(q). 

As usual in the DMRG approach, we consider open 
chains. The numerical calculations were performed u sing 
the dynamic block-state selection (DBSS) approach [29|. 
The threshold value of the quantum information loss x 
was set to 10~ 8 for the spin models and to 10 -4 for the 
fermionic model, and the upper cutoff on the number of 
block state was set to Af max = 1500. 

At and below the VBS point, i.e., for -tt/4 < 8 < 
#vbSj s (0 increases with I for small I, saturates due to 
the Haldane gap [3(J > an d then goes down to zero again as 
I approaches N. The Fourier spectrum s(q) is a smooth 
function of q (except for the q = component). The 
transformed entropy s(q) at the VBS point, depicted in 
Fig. illustrates this behavior. For 8 slightly larger than 
$vbs, however, we find that s(l) does not increase to 
the saturation value purely monotonically. Instead, an 
incommensurate oscillation is superimposed. For some- 
what larger 8 values, 8 > 0.137T, this oscillation persists 
in the saturated region, i.e., for blocks much longer than 
the correlation length. This leads to a new peak in s(q) 
which moves from small q towards q = 2n/3 as the ULS 
point is approached, and gets larger and narrower, as can 
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be seen in Fig. This 9 value is slightly smaller than, 
but close to, the Lifshitz point. 
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FIG. 3: Fourier-transformed entropy s(q) for various (9, ob- 
tained for a finite chain with N = 180 lattice sites. The lines 
are guides to the eyes. The 5 = point, which has a large 
positive value, is not shown. 

The spin-1/2 frustrated Heisenberg chain, Hamiltonian 
{J2J, is also known to develop incommensurate correla- 
tions, for values of J'/ J > 0.5 (the MG point). As 
shown in Fig.0] the entropies of blocks of length N/2 and 
N/2 + 1, although substantially different in value, both 
display a minimum as a function of J'/ J at the MG point. 
Thus, the transition from commensurate (dimerized) to 




FIG. 4: Entropy of blocks of length N/2 and N/2 + 1 as a 
function of J' /J for the frustrated spin-1/2 Heisenberg chain 
for various chain lengths. 

incommensurate correlations again is marked by an ex- 
ternum of the block entropy. Similarly to the behavior 
in the neighborhood of the VBS point, in this model a 
long wavelength oscillation appears in s(l) right above 
J' / J > 0.5, which leads to a new peak in the Fourier 
spectrum. The displacement with increasing J'/J is in 
this case, however, towards q — tt/2, as can be seen in 
Fig. [S] This peak becomes larger and narrower for in- 
creasing J 1 / J , and again is a signature of incommensu- 
rability. 

It is also interesting to examine the behavior of the 
block entropy in the lowest-lying triplet excited states. 
This is shown in Fourier-transformed representation, 
s(q), in Fig. for several J'/J values. At J' = 0.5 J, the 
peak at q = ir is displaced to q = 7f(l — 1/N). As J'/J 
increases, we find two oppositely moving peaks. One ap- 
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FIG. 5: Fourier-transformed entropy s(q) for various values of 
J' J for the frustrated Heisenberg chain with N — 128 lattice 
sites. 



pears exactly where the peak was found for the ground 
state, while the other occurs at ir — q. By also calculating 
the structure factor, S(q), we found that this second peak 
is located at the same (J'/J-dependent) wave vector at 
which S(Q) has its maximum |3l| . For Hamiltonian J2J, 
the block entropy of the first triplet state also has two 
oppositely moving peaks. One is again at the same lo- 
cation as the incommensurate peak in the ground state, 
while the other is at n — q/2. Both peaks move towards 
2n/3 with increasing 8. 
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FIG. 6: Fourier-transformed entropy 3(g) for various values of 
J'/J for the triplet excited state of the frustrated Heisenberg 
chain. The arrows indicate the direction of movement of the 
peaks with increasing J'/J. 

Having demonstrated the usefulness of studying the 
entropy profiles for models where the quantum critical 
points are known, we now use this procedure to study 
an commensurate-incommensurate transition in the ID 
t — t' — U Hubbard model 



(5) 



^ _ ( C lcT C 'i + l(T + C l+l(T C i, 
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which has been investigated recently [3^. For the half- 
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filled case (and setting t = 1), the competition between 
if and the Coulomb energy U will determine whether the 
system is an insulator (t' < t' s ) or a metal [t' > t' c ). For 
finite U values, the transition between these two states 
occurs in two steps. First the spin gap opens at t' s and 
then the charge gap closes at a larger value t' c . Between 
these two points, the wave vector becomes incommen- 
surate for t' > t' lc — the commensurate-incommensurate 
transition is independent of the metal-insulator transi- 
tion 0. 

As expected for a commensurate-incommensurate 
transition, we find that the entropy of blocks of length 
N/2 and N/2+1 display an extremum as a function of t'. 
For very large U values, where the model is equivalent to 
the frustrated spin-1/2 Heisenberg chain, the extremum 
occurs at t' IC ~ l/\/2, which maps to the MG point. For 
t' > tj C , an incommensurate oscillation in s(l) becomes 
apparent. This behavior can be seen in Fig. [7| for {7 = 3, 
a value chosen so that our results can be directly com- 
pared to those of Rcf. 32] . When s(q) is analyzed, it 
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FIG. 7: Block entropy profile of a finite chain with N = 64 
sites for the t — t' — U Hubbard model for U = 3 at t' = 
0.6, 0.625, 0.65, 0.675 (from top to bottom). 

is found that a new peak appears in the spectrum and 
again moves from small q towards q = tt/2 with the am- 
plitude of s(tt) decreasing with increasing i! . Therefore, 
the commensurate-incommensurate phase boundary can 
be easily determined by finding the extrema of s(N/2) 
as a function of i' for various U values. The detailed en- 
tropy analysis of the full two-dimensional phase diagram 
will be presented in subsequent work. 

In conclusion, we have shown that the length depen- 
dence of the block entropy and its Fourier spectrum, de- 
termined for finite systems, can be used to characterize 
phases in which the correlation function has an oscilla- 
tory character. In addition, an extremum in the block 
entropy as a function of the relevant model parameter, 
which, in general, signals the appearance of or change in 



a symmetry in the wave function, can also correspond 
to disordered points. In this case, however, the entropy 
curve does not show anomalous behavior because this is 
not a phase transition in the conventional sense. When 
the decaying correlation function has an incommensu- 
rate oscillation, a new peak appeares close to q = in 
the Fourier spectrum and moves towards a commensurate 
wave vector as the control parameter is adjusted. In the 
entropy of the triplet excited states, in addition to a peak 
at the same position as in the ground-state entropy, an- 
other peak appears at the wave vector of the peak in the 
static structure factor. Our method is ideal for use in con- 
junction with the density-matrix renormalization-group 
algorithm because the block entropy profile is generated 
as a by-product of the DMRG procedure. This allows 
the more difficult and sensitive calculation of correlation 
functions to be avoided. 
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